Scaling of the F2 structure function in nuclei and quark distributions at x > 1 
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We present new data on electron scattering from a range of nuclei taken in Hall C at Jefferson Lab. 
For heavy nuclei, we observe a rapid falloff in the cross section for x > 1, which is sensitive to short 
range contributions to the nuclear wave-function, and in deep inelastic scattering corresponds to 
probing extremely high momentum quarks. This result agrees with higher energy muon scattering 
measurements, but is in sharp contrast to neutrino scattering measurements which suggested a 
dramatic enhancement in the distribution of the 'super-fast' quarks probed at x > 1. The falloff at 
a: > 1 is noticeably stronger in and ''He, but nearly identical for all heavier nuclei. 

PACS numbers: 13.60.Hb, 24.85.-(-p, 25.30.Fj 



The quark structure of nuclei is extremely complex, 
and a detailed understanding of nuclei at the quark level 
requires the careful incorporation of nucleonic degrees 
of freedom and interactions as well as the dynamics of 
quarks and gluons. In inclusive electron scattering from 
nuclei, the cross section is characterized by the structure 
functions Fi{v, Q^) and F2{v, Q^), where i/ is the energy 
transfer and —Q^ is the square of the four-momentum 
transfer. At high Q^, the reaction is dominated by 
elastic scattering from quasi-free quarks, and one can 
probe the momentum distribution of the quarks. In 
the Bjorken limit (z/, — > oo), the quark mass and 
transverse momenta are negligible compared to the en- 
ergy and momentum of the probe, and the scattering 
is sensitive only to the the quark longitudinal momen- 



tum, where x — / {2Mv) is the fraction of the hadrons 
longitudinal momentum carried by the quark in the infi- 
nite momentum frame. In this deep inelastic scattering 
(DIS) limit, the structure functions exhibit scaling, i.e. 

ming independent of at fixed 
X, with F2{x) proportional to a charge- weighted sum of 
the quark momentum distributions in the target. 

As one moves away from the Bjorken limit, there are 
deviations from perfect scaling. At finite-Q^, kinematical 
corrections yield a dependence that can be large for 
low or large x. While these scaling violations have his- 
torically been called "target mass" corrections [H, they 
are in fact independent of the mass of the target for a 
quark of a given longitudinal momentum. At lower , 
there are also significant contributions from higher twist 
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effects, e.g. structure due to quark-quark and quark- 
gluon correlations which appear most clearly as strong 
resonance structure. These scaling violating terms make 
extraction of the quark distributions most straightfor- 
ward at high energies. QCD evolution yields an approx- 
imately logarithmic dependence at all values, but 
this is a true scale dependence of the parton distributions. 

The early expectation was that the nuclear quark mo- 
mentum distribution would be a convolution of the distri- 
bution of nucleons in a nucleus with the the distribution 
of quarks in the nucleons. Contrary to these expecta- 
tions, measurements of inclusive scattering from nuclei 
showed a 10-20% suppression of high momentum quarks 
(0.3 < a; < 0.8) in heavy nuclei demonstrating that 
the quark distributions in nuclei are not simply a sum of 
the proton and neutron's quark distributions. 

The quark distribution at a; > 1 is extremely small in 
the convolution model, as the nucleon quark distributions 
fall rapidly as x 1 and there are very few fast nucle- 
ons available to boost the quarks to a: > 1. The bulk of 
these 'super- fast' quarks come from nucleons above the 
Fermi momentum, which are generated by the strongly 
repulsive core of the N-N interaction; they reflect the 
short-range correlations (SRCs) in the ground state nu- 
clear wave- function 0,0]. Exotic configurations, such as 
6-quark bags, may provide an even more efficient mecha- 
nism for generating very high momentum quarks, as the 
quarks from two nucleons can more freely share momen- 
tum [1, 0] . It is clear that a holistic explanation of DIS 
from nuclei must describe the behavior of the structure 
functions in the full kinematic range, and measurements 
at a; > 1 are necessary to illuminate the presence of short 
range structure in nuclei. 

There are only two high energy measurements thus 
far for a; ^ 1, and they yield dramatically different re- 
sults. Muon scattering data from the BCDMS collabora- 
tion [H for Q2 from 52-200 GeV^ show a rapid falloff in 
^2(3;). They find ^2(2^) exp(— s • x) with s = 16.5 ± 0.6 
for 0.75 < a: < 1.05, which suggests the need for rel- 
atively modest contributions from short-range correla- 
tions. Neutrino scattering measurements from the CCFR 
collaboration at = 125 QeV^ yield s = 8.3 ± 0.7 
for 0.75 < a: < 1.2, which has been interpreted as an 
indication of exceptionally large strength from short- 
range correlations or the need for other, more exotic, 
contributions. However, both measurements have impor- 
tant limitations: CCFR had to make significant correc- 
tions due to the poor resolution on their reconstructed x 
value, while BCDMS was only able to extract Fi up to 
X = 1.05. It is unclear if this is sufficient to make mean- 
ingful comparisons to model predictions of short range 
structure in nuclei, as this is not expected to dominate 
until X ^ 1.2 ^, |8||- More extensive measurements have 
been made at lower energies, but they have been limited 
to a; « 1 1, [3 or Q2 < 5 QeV^ ^M^- 

We present the results of JLab E02-019, which made 



new measurements of inclusive scattering from nuclei, 
covering an expanded range in x and Q^. Data were 
taken for few-body and heavy nuclei, covering both the 
region of the EMC effect [Tl and X > 1. The measure- 
ment was performed in Hall C at the Thomas Jefferson 
National Accelerator Facility in 2004. A continuous wave 
electron beam of 5.766 GeV and current of w80 /iA was 
supplied. Electrons scattered from the target were de- 
tected using the High Momentum Spectrometer (HMS) 
at 61 = 18°, 22°, 26°, 32°, 40°, and 50°, corresponding 
to 2 ^ ^ 9 GeV^. Data were taken on cryogenic ^H, 
■^He, *He targets, sohd Be, C, Cu, and An targets, as 
well as Al targets used to measure and correct for the 
contribution from the walls of the cryogenic target cells. 

Electrons were selected using the HMS gas Cerenkov 
and electromagnetic calorimeter detectors with efficien- 
cies of >98% and >99.7%, respectively, and negligible 
pion contamination. Data were also taken with the HMS 
in positive polarity to determine the contribution from 
charge-symmetric processes. The systematic error asso- 
ciated with positron subtraction was negligible except at 
the largest values, where it is under 2%. The uncer- 
tainty due to the spectrometer acceptance is 1.4%. 

The largest sources of systematic uncertainty at high 
X come from the beam energy (0.05%), HMS central mo- 
mentum (0.05%) and angle (0.5mr) settings. The impact 
of these on the cross section is typically small (1-2% for 
a: < 1) for all but the largest x values, where the uncer- 
tainty can reach 4-6%. The systematic uncertainty that 
arises from subtraction of the aluminum end-caps (cryo- 
genic targets only) also grows with increasing x values, as 
does the relative contribution of the aluminum end-caps 
to the measured cross section, giving a range in error of 
0.3-2.4%. The cross sections also had to be corrected 
for bin-centering, radiative, and Coulomb effects. The 
systematic errors associated with those corrections are 
0.5%, 1.4%, and <2%, respectively. Details on the anal- 
ysis and uncertainties, as well as the various corrections 
applied can be found in Ref. 15|. The total systematic 



uncertainty on the extracted cross sections is below 4% 
for a; < 1, and up to 6% in the 1 < a: ^ 2 range. 

The structure function per nucleon, F2{v,Q^), is ex- 
tracted from the measured cross section as follows: 



d^'^E a™oti[l + 2 tan2(0/2) 



l + R 



(1) 



with i?=(TL/crT=(0.32 GeV2)/Q2 [9] and 5R/R^bQ%, 
yielding an additional uncertainty of ;$1% in i^2- 

Scaling of the proton structure functions has been ob- 
served over a large kinematic range in high energy inclu- 
sive scattering. Based on these measurements, the DIS 
region is often taken to be W'^ > 4 GeV^, > 1 GeV^, 
where — Mp + iMpU — is the square of the 
invariant mass of the undetected hadronic system and 
Mp is the proton mass. It has been shown that scaling 
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violations are reduced when one examines F2 at fixed 
^ = 2x/ (1 + r) [1^, where r = ^^1 + jv"^. ^ is equiva- 
lent to X in the Bjorken limit, but when examining scal- 
ing at fixed ^, rather than fixed x, the observed scaling 



behavior extends to lower W"^ [12|, |13| , corresponding to 
larger ^ values. This improved scaling can be seen clearly 
in Fig. [1] where the upper sets of curves show F2 for car- 
bon plotted against x (red squares) and f (green circles) 
over a range of values. The extended ^-scaling of the 
nuclear structure function, seen to begin above 3-4 GeV^ , 



may allow access to quark distributions for ^ ^ 1 [ITlll 
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FIG. 1: (Color online) F2 for the E02-019 carbon data (2 ^ 
^ 9 GeV'^) as a function of x (top) and ^ (middle), and 
Fj"^ vs ^ (bottom). In each case, the higher points correspond 
to the smaller scattering angles (lowest values). 

We compare our data to higher measurements, us- 
ing a partonic framework to look for deviations from the 
scaling picture. Rather than simply examining F2 as a 
function of ^, as was done previously, we account for the 
kinematical scaling violations using the prescription of 
Ref. fl| (Eq.(23)) and study the scaling of F^°\^,Q^y. 

F^"\i,Q') = F2{x,Q')-^^h2{^,Q') 



Q 

(0), 



4„5 



52(^,< 



(2) 



and we obtain a fit for F2^\^, Qq) from a subset of these 
data, chosen to minimize contributions from quasielastic 
scattering. This simple fit provides a reasonable descrip- 
tion of the global data set (see Fig. [5]), with deviations at 
low Q^, in particular near the quasielastic peak (^ « 0.85) 
and for the largest values of The /12 and g2 terms yield 
corrections of up to 15% at the lower Q-^ values, but ^5% 
for > 5 GeV^. We estimate the model dependence in 
the extraction of Fo""* to be <,2%. 
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FIG. 2: (Color online) F^"'' vs for fixed ^ value. For 
this work and BCDMS, the carbon data are shown, while 
the SLAC points are carbon pseudo-data taken from mea- 
surements on deuterium. The solid curves are the global fit, 
the short horizontal red lines show the BCDMS ^=1.15 upper 
limit, and the green crosses show the falloff between ^=0.75 
and ^=1.05 based on the CCFR data (see text for details). 



where /i2(e,Q') = u-^F^''>{u,Q^)du and 52(e,Q') = 
jf' v-'^{v~0F2°\'"^Q'^)d'"- One could incorporate these Figure [T] shows F^"^ vs £_ (blue triangles), which has 



effects into a partonic model for F2 , rather than extract- 
ing an "idealized" scaling function, but this approach 
minimizes the dependence in the scaling function, 
making it easier to directly compare different data sets. 

To calculate /12 and 52, we use a factorized model for 
-^2"' Q^)' with a common dependence for all targets 
and a simple fit to F2°\^,Qq) for each nucleus. In the 
partonic picture, the dependence should come only 
from QGD evolution. Because we cannot determine QCD 
evolution without knowing the partonic structure, we fit 
the the dependence of the world's data (shown in 
Fig. [2]), excluding our lower Q"^ points, at several ^ values 
to a functional form chosen to be consistent with evolu- 
tion. This is used to scale our new data to Qq = 7 GeV^ , 



somewhat greater dependence at large ^ than the un- 
corrected structure function F2 (green circles). The main 
difference between ^-scaling and our scaling analysis is 
the factor x'^ /{(^"^r^) in front of the leading term in Eq. [21 
as the /i2 and 32 terms are relatively small. Neglecting 
this pre-factor introduces an additional dependence 
that approximately cancels that of the QCD evolution, 
resulting in an artificially small dependence in the 
naive ^-scaling analysis. 

Figure [2] shows the carbon results for F2^\^,Q^), 
scaled to fixed values of ^ using our global fit. The SLAC 
points are deuterium data [19|, multiplied by the SLAC 
E139 ^20] fit to the carbon-to-deuteron structure function 
ratio, yielding carbon pseudo-data to provide a continu- 
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ous range for lower ^ values. For all data sets, F^^ 
is extracted from the measured structure functions using 
the global fit to calculate 52 and /i2- For ^ < 0.75, where 
the high data determine the evolution, our data are 
in excellent agreement with this dependence down to 
= 3 GeV^. The observed dependence grows slowly 
with ^ over this region, and with a continued increase at 
higher ^ values, our highest measurements are con- 
sistent with SLAG and BCDMS. For large ^ values at 
low Q^, our data deviate from this dependence due 
to higher twist contributions, especially in the vicinity of 
the quasielastic peak k, 0.85) 

The CCFR measurement did not explicitly extract val- 
ues of F2, but obtained a fit to the falloff at large ^. We 
illustrate this falloff by normalizing to our global fit at 
^ = 0.75 and applying the CCFR ^ dependence to extract 
plf^ at ^ = 0.75, 0.85, 0.95, and 1.05, shown as green 
crosses. This behavior is clearly inconsistent with the 
overall behavior of the structure function extracted from 
electron and muon scattering. The BCDMS data exhibit 
somewhat unusual behavior at large ^. Above ^ = 0.65, 
the BCDMS data show little or no dependence, even 
though one expects noticeable QCD evolution. 




FIG. 3: (Color online) The slope s in exp(— s-^) as a function 
ofQ^. The targets are offset in for visibility. Open symbols 
for and ^He at low are cases where the kinematic limit 
for the nucleus [x ~ A) corresponds to ^ ^ 1.25. 

To quantitatively examine the falloff of our structure 
function at large ^, we perform a fit similar to BCDMS 
and CCFR. We take the data from a fixed scattering 
angle, use the global fit to interpolate to a fixed value 
(corresponding to ^ = 1.1), and fit ^2*^^ (^, Q^) to exp(— s- 
^) for 1 < ^ < 1.25. The lower ^ limit is chosen to avoid 
regions where the quasielastic peak leads to noticeable 
deviations from scaling at low Q^, and the upper ^ limit 
is chosen so that there are data covering the full ^ range 
for all targets and values. We take the slope extracted 
from the 40° data ((5^=7.35 GeV^) as the main result, as 
this is the largest value with high statistics over the 
full ^ range. Data at smaller angles are used to examine 



the dependence of the result. 



TABLE I: Extracted values of the slopes for all nuclei. The 
uncertainties includes statistics and systematics; the latter 
are typically ~0.4 and dominate the uncertainty. 



A 


= 2.79 


3.75 


4.68 


5.95 


7.35 


2 


14.7±0.6 


16.6±0.5 


17.1±0.5 


17.5±0.7 


16.8±0.8 


3 


15.1±0.6 


15.6±0.5 


16.6±0.6 


16.4±0.7 


17.4±0.9 


4 


14.8±0.5 


14.7±0.5 


15.1±0.5 


15.1±0.5 


14.8±0.6 


9 


14.7±0.5 


14.6±0.4 


14.8±0.4 


15.0±0.4 


15.0±0.5 


12 


14.7±0.5 


14.3±0.4 


14.9±0.4 


14.5±0.4 


15.0±0.5 


64 


13.8±0.5 


13.7±0.4 


14.2±0.4 


13.9±0.4 


14.1±0.5 


197 


13.5±0.6 


13.6±0.5 


13.9±0.5 


13.9±0.5 


14.3±0.7 



The extracted slopes are shown in Table I] and Fig. [3l 
Above 4 GeV^, there is no systematic dependence, 
and at lower Q^, only the and '^He results change 
significantly. We observe nearly identical behavior in the 
high-^ falloff for all nuclei except and ^He, which have 
a larger slope and thus a steeper falloff with ^. 

We obtain s = 15.0 ± 0.5 for carbon, s = 14.1 ± 0.5 
for copper (our closest nucleus to the CCFR iron tar- 
get), showing that the large difference between BCDMS 
and CCFR is not related to the difference in target nu- 
clei. Note that BCDMS and CCFR extract slopes from 
^2(2^) instead of F^\(^), although the difference would 
increase their slopes by less then 0.5 (0.1) for the BCDMS 
(CCFR). Further complicating direct comparison is the 
fact that none of these experiments cover the same ^ 
range. For our new data, variations in the ^ limits of 
0.05-0.1 can change the extracted slope by 0.5-1.0. 
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FIG. 4: (Color online) The extracted scaling structure func- 
tion per nucleon for all nuclei at (5^=2.79 and 7.35 GeV^. 
The deuteron data are kinematically limited to a; < 2, corre- 
sponding to ^ ^ 1.15 for the low setting. 
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We have focused on heavier nuclei for comparison to 
BCDMS and CCFR, but have also obtained a signifi- 
cantly expanded body of data for light nuclei, shown in 
Fig. 131 At low Q^, the light nuclei show a clear quasielas- 
tic peak, while at higher the peak is almost entirely 
washed out. For heavier nuclei, the extracted scaling 
function per nuclcon is nearly identical at all ^ values. 
However, for and "^He, there is a significant reduction 
in strength for ^ ^ 1, which is observed at all values. 

In summary, we have made extensive measurements 
of the large ^ structure functions for from 2- 
9 GeV^. We have extracted the scaling structure func- 
tion, -F2°^(C, Q^), and shown that it is consistent with 
a nearly logarithmic dependence over a significant 
range of ^ and Q^. These new data do not show a need 
for extremely large contributions from short-range cor- 
relations or other, more exotic, short range structures, 
as suggested by the CCFR result. The large ^ behavior 
of our data is consistent with the BCDMS results, but 
our results extend to significantly higher ^ values, where 
one expects to be most sensitive to short range struc- 
ture U Q . The new data, covering a range of nuclei at 
large ^, can be used to directly constrain calculations of 
the effect of short-range correlations or multi-quark con- 
figurations in our kinematic regime for a wide range of 
nuclei. A future 12 GeV JLab measurement 2l| will dou- 
ble the range for these ^ values, moving to a region 
where we can directly extract the parton distributions of 
the super-fast quarks in nuclei. 
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